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Current understanding of the brittleness of glass is limited by our poor understanding and control over the 
microscopic structure. In this study, we used a pressure quenching route to tune the structure of silica glass 
in a controllable manner, and observed a systematic increase in ductility in samples quenched under 
increasingly higher pressure. The brittle to ductile transition in densified silica glass can be attributed to the 
critical role of 5-fold Si coordination defects (bonded to 5 O neighbors) in facilitating shear deformation and 
in dissipating energy by converting back to the 4-fold coordination state during deformation. As an 
archetypal glass former and one of the most abundant minerals in the Earth's crest, a fundamental 
understanding of the microscopic structure underpinning the ductility of silica glass will not only pave the 
way toward rational design of strong glasses, but also advance our knowledge of the geological processes in 
the Earth's interior. 

Silica glass (a-sUica), owing to its low atomic packing density, can be permanently densified to a large extent 
(>20%) by static cold compression at room temperature''^, shock compression' or by pressure quenching 
from melt*. Extensive experimental and computational studies have been conducted to understand the 
structure and properties of o-silica under static cold compression at room temperature^ '". Raman spectra of 
shock-retrieved o-silica (experienced peak pressures of 52 GPa and temperatures of 3000 K) show substantial 
differences from those of a statically cold compressed sample'. This indicates that the compression mechanisms 
in a-silica at low temperature seem to be different from those at high temperature, which may be closely related to 
the compression mechanisms in the melt. This view was supported by NMR studies in silicate glasses quenched 
from melts under high pressure, in which appreciable amount of 5- and 6-fold coordinated Si was identified, while 
only 4-fold Si atoms were observed in normally melt-quenched samples^""^'. Our previous study showed that 
elastic properties of o-silica can be tuned by the pressure quenching route in a controllable manner''. It is the goal 
of this work to systematically study the plastic deformation mechanisms of pressure-quenched a-silica. Direct 
probing the microscopic changes associated with the plastic deformation in glass is a formidable challenge for 
experiment''. Here we resorted to large-scale molecular dynamics (MD) simulation in this study, which has an 
advantage over experiment in that it can provide atomistic-level structural details associated with the plastic 
deformation during mechanical tests. 

A modified BKS potentiaP'''^^ implemented in the LAMMPS^*" package was used to prepare a series of pressure- 
quenched a-sUica samples (details can be found in the Methods part and in Supporting Information (SI)). 
Densified a-silica samples are labeled by the pressure applied during the quenching process for convenience. 
In this study, we investigated the deformation mechanisms of densified a-silica using three large-scale mechanical 
tests (i.e., uniaxial tension test, V-crack tension test and nanoindentation test) at the same homologous temper- 
ature of 0.1 Tf (fictive temperature, defined and calculated in SI). The structural origins responsible for the 
enhanced ductility in densified a-sUica were attributed to 5-fold Si coordination defects formed during the 
pressure quenching process. The amount of 5-fold Si coordination defects can be controlled by the quenching 
pressure, which in turn determines the deformation modes of densified a-sUica. As an archetypal glass-former 
and one of the most abundant minerals in the Earth's crust, a fundamental understanding of the microscopic 
structure underpinning the ductility of a-sUica wUl not only pave the way toward the rational design of strong 
glasses''' but also shed light on the geochemical and geophysical processes in the Earth's interior' "*"" ''^. 

Results 

A range of constant hydrostatic pressures were applied during the melt-quenching process to prepare a series of 
densified samples'*. Upon pressure release at 0.1 Tf, Fig. 1(a) shows that the fraction of 5-fold Si atoms steadily 
increases with the quenching pressure (more than 20% in the 15 GPa sample), which constitute the major 
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Figure 1 | Coordination defects in pressure-quenched a-silica. (a) Fraction of 4-, 5- and 6-fold Si versus quenching pressure; (b) Configuration of 
two types of Si atoms (red: 4-fold, green: 5-fold) in the 15 GPa sample, O atoms are omitted for clarity; (c) Si-Si radial distribution function of 4-fold and 
5-fold Si atoms in the 0 GPa and 15 GPa sample (represented by the subscript 4 and 5 in g(r), respectively). 



coordination defects in densified o-silica. There is only ~ 1 % 6-fold Si 
coordination defects in the 15 GPa sample, which indicates that sUica 
melt/glass seems to compress differently from the crystalline coun- 
terparts under pressure where the coordination number of Si changes 
from 4 to 6^*'. Distribution of the 4- and 5 -fold Si in the 15 GPa sample 
is shown in Fig. 1(b), indicating a tendency of cluster formation in 5- 
fold Si atoms. Radial distribution functions of the 4-fold and 5-fold Si 
in the 0 GPa and 15 GPa sample in Fig. 1(c) confirm that some of the 
5-fold Si atoms tend to stay closer together, indicated by the inner 
g(r)55 peak below 3 A. Clusters of 5-fold Si atoms are uniformly 
distributed throughout the sample, which may correspond to the 
high density regions in the original silica melt". 

To investigate the deformation mechanisms of densified a-sUica, 
we first carried out uniaxial tension tests in as-quenched samples (0, 



2, 4, 6, 8, 10, 12 and 15 GPa). To avoid the system size effect as 
reported in our previous study'", dimensions of samples for uniaxial 
tension tests have been properly chosen (see Table. S2). In Fig. 2(a), a 
clear brittle to ductile transition can be found in samples quenched 
under pressures of 2 to 4 GPa. The 0 GPa sample (i.e., pristine 
a-sUica) displays a clear brittle fracture behavior, consistent with 
observations from previous experiments-"'*" and simulations"*''"'''^. 
Above 4 GPa, all densified samples don't fracture even up to 50% 
uniaxial strain, indicating a typical ductile behavior. A local shear 
strain analysis'*'' was performed at selected tensile states during uni- 
axial tension tests for the 0 GPa and 8 GPa sample (see Fig. 2(b)). In 
general, a higher local shear strain value (i.e., more reddish) means 
that a local cluster has been heavily shear deformed. For the 0 GPa 
sample, void formation and crack propagation mechanism can be 
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envisaged from the local shear strain mapping in Fig. 2(b). Voids 
probably form in heavily sheared regions due to decohesion resulted 
from local heating"" or due to heterogeneous nanoscale cavitations'*^. 
In contrast, for the 8 GPa sample, despite a large portion of the 
sample is heavily shear deformed, no such voids form, implying that 
void formation has a higher activation energy than shear deforma- 
tion in densified a-sUica. 

To study the resistance of densified a-silica to propagation of pre- 
existing cracks, we carried out a series of V-crack tension tests in the 
0, 8, 10, 12 and 15 GPa samples. Fig. 3 shows that with increasing 
quenching pressure, samples display more and more plastic 
deformation before fracture, evidenced by shear banding and crack 
blunting. The strain at fracture or fracture strain, as an indicator of 
ductility, increases from 10.5% to 24.1% when the quenching pres- 
sure increases from 0 to 15 GPa (see Fig. S13). For all samples, pre- 
existing cracks always induce stress concentration around crack tip 
regions and can serve as nucleation sites for new voids. However, for 
the 0 GPa sample, the crack never blunts, and the void formation and 
coalesce phenomena at crack tip, although regarded as indication of 
nanoscale ductility*'', would not lead to a macroscopic ductile beha- 
vior due to high crack propagation speed and localized void forma- 
tion very close to the crack tip'*" '". For densified samples, the crack 
substantially blunts before propagating together with extensive shear 



flow events (as evidenced by the local shear strain mapping in the 15 
GPa sample in Fig. 3(b)). As seen in Fig. S12, the fracture surface is 
very smooth for the 0 GPa sample, similar to Guin and Wiederhorn's 
AFM results'*'; it becomes increasingly rougher with the increase of 
quenching pressure. The fracture surface roughness is 7.6, 20.2 and 
33.3 A for the 0, 8 and 15 GPa sample, respectively, another signature 
of enhanced ductility in densified samples. 

In addition, nanoindentation tests in the 0 GPa and 15 GPa sample 
were performed using a V-shaped nanoindenter with a round tip. 
Load- displacement curves of the two nanoindentation tests are plot- 
ted in Fig. 4(a) and configuration snapshots after a fuU unloading of 
the indenter are compared in Fig. 4(b). As seen in Fig. 4(a), the 
densified a-silica (i.e., the 15 GPa sample) definitely has a higher 
contact elastic modulus and nano-hardness than the pristine a-sUica 
(i.e., the 0 GPa sample). The serrated stress drops in the plastic region 
(after —1.0 nm indentation depth) can be attributed to permanent 
local rearrangements (e.g., shear flow events'*"*). By comparing the 
after-unloading snapshots, the densified a-silica shows much 
more pile-up than the pristine one in Fig. 4(b). It can be seen in 
Fig. 4(b) that the densified a-silica tends to deform predominantly 
through shear flow instead of permanent densification under sharp 
contact loading. The shear flow under nanoindentation leads to more 
pile-up after unloading, in excellent agreement with Rouxel et al. 's 
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Figure 2 | Deformation of pressure-quenched a-silica during uniaxial tension test, (a) Stress-strain curves ofuniaxial tension test for pressure-quenched 
samples (0, 2, 4, 6, 8 and 10 GPa); (b) Local shear strain mapping of uniaxial tension test for the 0 GPa (top) and 8 GPa (bottom) sample. 
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Figure 3 | Deformation of pressure-quenched a-silica during V-crack tension test, (a) Stress-strain curves ofV-crack tension test for pressure-quenched 
samples (0, 8, 10, 12 and 15 GPa); (b) Local shear strain mapping of V-crack tension test for the 0 GPa (upper), 8 GPa (middle) and 15 GPa (bottom) 
sample. 



experimental observations in densified a-silica obtained from static 
cold compression at room temperature^. While in the pristine a- 
silica, permanent densification is the dominant deformation mode 
under indentation due to its low atomic packing density (Fig. 4(b)), 
consistent with Yoshida's experimental observations'*''. 

Discussion 

A simple model based on the competition between two fracture 
mechanisms, namely "shear banding" versus "cleavage", is used 



(more details in SI) to understand the brittle to ductile transition 
in the densified a-silica from the energy point of view. The fracture 
energy for the cleavage mode is 2Yfse X A and for the shear banding 
mode is Ysbe / cos(6) X A, where A is the cross section area perpen- 
dicular to the tensile axis, Yf^e is fracture surface energy per unit area, 
9 is the shear band angle and taken as 45° for simplicity in this study, 
and Ysbc is shear banding energy per unit area (assuming the thick- 
ness of shear bands in uniaxial tension test and pure shear test are the 
same, see SI). Since the cross section A is the same for two different 



SCIENTIFIC REPORTS | 4:5035 | DOI: 1 0. 1 038/srep05035 



4 




0 5 10 15 

Displacement of indenter (nm) 



OGPa 15GPa 



Atomic configuration 



Local shear strain 



Local density 



Local density change 




Figure 4 | Deformation of pressure-quenched a-sUica during nanoindentation test, (a) Load-displacement curves of nanoindentation test for tlie 0 GPa 
and 15 GPa sample; (b) Atomic configuration (first row), local shear strain (second row), local density (third row) and local density change 
mapping under indenter (fourth row) after nanoindentation test for the 0 GPa and 15 GPa sample. 



fracture modes, the net comparison is between 2yfse and ysbe / cos(9). 
If 2Yfse > ysbe/cos(6), shear banding fracture will dominate; other- 
wise, cleavage mode will be more feasible. The crossover of the two 
energies in Fig. 5(a) around 2-4 GPa implies the fracture mode 
changes from cleavage to shear flow, consistent with the trend 
observed in Fig. 2(a). 

In MD simulations, we can directly observe the dynamic evolution 
of the 5-fold coordination defects during large-scale mechanical tests. 
The fraction of the 4- and 5 -fold Si atoms in uniaxial tension tests are 
plotted for the 0, 4, 8 and 15 GPa samples in Fig. 5(b) and very 
different behaviors can observed in samples with different amount 
of 5-fold Si coordination defects after the pressure quenching process. 
For the 0 GPa sample with negligible amount of 5-fold Si, it exhibits a 
clear brittle fracture at the maximimi stress. For samples densified 
above 8 GPa, the 5-fold Si converts back to the 4-fold coordination 
state upon initial uniaxial tension, which gives the steady decrease in 
the first region in Fig. 5(b), then approaching a constant population 
of —4% at larger strains. In the 4 GPa sample, the initial population of 
the 5-fold Si coordination defects is about 4%, which decreases 
slightly upon uniaxial tension then approaches the steady state as 



others. It is very interesting to see in Fig. 5(b) that as long as a sample 
has enough 5 -fold coordination defects to start with, it can reach a 
steady flow state as shown in Fig. 2(b). This converging behavior is 
very similar to what was observed in simple shear tests of model 
metallic glasses^". The easiness of shear flow in densified a-sUica as 
evidenced in aU three mechanical tests can be attributed to a less rigid 
local environment associated with 5-fold Si coordination defects, 
more amenable to local shear deformation as compared with the 4- 
fold coordination state as shown in Fig. S20-S21. 

Fig. 5(b) also shows that the coordination state conversion in 
densified a-silica provides a new deformation mode besides the bond 
stretching upon uniaxial tension. Fig. S23-25 further demonstrate 
that the conversion from the 5-fold (dense packing) to the 4-fold 
(loose packing) state provides additional mechanism for energy dis- 
sipation in densified fl-silica. To some extent, the two-level deforma- 
tion mechanism in densified a-silica is similar to that in polymers: 
untwining the network chains before bonds in the network are 
stretched. In the pristine a-sUica, the dominant energy dissipation 
mechanism upon expansion is simply to stretch bonds, leading to 
earlier bond rupture and cavitation. 
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Figure 5 | Origins of the enhanced ductility in densified a-silica. (a) Fracture surface energy and shear banding energy of fl-silica versus quenching 
pressure; (b) Fraction of 4- and 5-fold Si as a function of strain during uniaxial tension test for densified samples (0, 4, 8 and 15 GPa). (c) Snapshots of the 
8 GPa sample at different times at 35% strain in the uniaxial tension test. Note: red for 4-fold Si, green for 5-fold Si, O atoms are omitted for clarity. Inside 
the dashed circle, two 5-fold Si atoms appear then disappear, indicating the migration of 5-fold Si atoms inside the sample while maintaining a constant 
population as seen in (b). 
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5-fold Si coordination defects have been shown to be responsible 
for the diffusivity maximum in silicate melts^', and proposed as one 
of the possible transition states for the high-temperature growth of 
quartz from glass""^. Previous cold-compression study of a-silica in 
MD simulation by Liang et al. also showed that 5-fold Si coordination 
defects enhance local rebonding and relaxation, and therefore plas- 
ticity (densification) beyond the elastic limit^'. Liang et al. further 
proposed that 5-fold Si-activated diffusion is responsible for the 
minimum yield strength around 10 GPa observed in mechanical 
strength measurements in fl-sUica^'. Fig. 5(c) shows that 5-fold Si 
atoms are migrating but maintaining a constant population during 
the flow state of the uniaxial tension test in Fig. 2(b). This indicates 
that 5-fold Si atoms act like "catalyst" for shear flow events or "plas- 
ticity carrier" analogous to the liquid-like component in amorphous 
Si proposed by Demkowicz and Argon""*, therefore playing an indis- 
pensable role in the enhanced ductility in densified a-silica. The 
negligible amount of 6-fold coordination defects observed in our 
pressure-quenched samples indicates that octahedral units such as 
those proposed by Stolper and Ahrens^ do not play an important role 
in the enhanced ductility in densified fl-sUica. 

In conclusion, densified a-silica prepared by the pressure-quench- 
ing route displays enhanced ductility evidenced by three large-scale 
mechanical tests: uniaxial tension test, V-crack tension test and 
nanoindentation test. Intrinsic fracture reflected by uniaxial tension 
test shifts from a brittle to a ductile mode as the quenching pressure 
increases and a simple model based on the cleavage and the shear 
banding fracture successfully explains the crossover in samples 
quenched under 2 — 4 GPa. Ductility of densified a-silica in the 
presence of V-shaped crack shows significant increase indicated by 
extensive shear banding and obvious crack blunting. There is much 
more pile-up under nanoindentation in densified fl-silica, indicating 
an easy activation of shear flow under indentation. Structural origins 
of the enhanced ductility in densified a-silica can be attributed to 5- 
fold Si coordination defects, which serve as "plasticity carrier" to 
promote shear flow and provide additional means for energy dissipa- 
tion by converting back to the 4-fold coordination state. Our study 
shows that the pressure quenching route may provide a novel way to 
quench the 5 -fold Si coordination defects, giving rise to enhanced 
ductility in densified a-silica. 

Methods 

Classical molecular dynamics (MD) simulations were conducted in the LAJVIMPS 
package''*' {http://lammps.sandia.gov/) using a modified version of the BKS poten- 
tiap* with a short range cutoff of 0.55 nm and a long-range Columbic cutoff of 1.0 nm. 
Previous studies have shown that this modified BKS potential can reproduce both 
bulk and surface properties of silica glass^^'"***. The Columbic interaction was calcu- 
lated via the Ewald summation technique with a relative precision of 10"^ in force. 
The Velocity- Verlet algorithm with a time step of 1.6 fs was used for integrating the 
equations of motion and the Nose-Hoover thermostat^^'^'' and barostat^^ were used to 
control the system temperature and pressure when necessary. 

A constant hydrostatic pressure was applied during the melt-quenching process 
(cooling rate of 10 K/ps), providing an additional driving force to density the glass 
sample. The pressure ramping and releasing rates during the pressure quenching 
process were fixed at 0.125 GPa/ps. The dimensions of pressure-quenched samples 
after pressure-release are listed in Table. SI and used in uniaxial tension tests. For V- 
crack tension tests and nanoindentation tests, they were duplicated in the x and y axis 
for multiple times. Periodic boundary conditions (PBCs) were applied in all direc- 
tions unless stated otherwise (e.g., in nanoindentation test). All mechanical tests were 
carried out at the same homologous temperature 0. 1 Tf (fictive temperature) to get rid 
of the temperature effect on the brittle to ductile transition observed in our previous 
study". 
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